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Abstract — We show how distributed control of reactive power 
can serve to regulate voltage and minimize resistive losses 
in a distribution circuit that includes a significant level of 
photovoltaic (PV) generation. To demonstrate the technique, 
we consider a radial distribution circuit with a single branch 
consisting of sequentially-arranged residential-scale loads that 
consume both real and reactive power. In parallel, some loads 
also have PV generation capability. We postulate that the 
inverters associated with each PV system are also capable 
of limited reactive power generation or consumption, and we 
seek to find the optimal dispatch of each inverter's reactive 
power to both maintain the voltage within an acceptable range 
and minimize the resistive losses over the entire circuit. We 
assume the complex impedance of the distribution circuit links 
and the instantaneous load and PV generation at each load 
are known. We compare the results of the optimal dispatch 
with a suboptimal local scheme that does not require any 
communication. On our model distribution circuit, we illustrate 
the feasibility of high levels of PV penetration and a significant 
(20% or higher) reduction in losses. 
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I. Introduction 

Utilities employ various equipment to control the oper- 
ation of primary distribution systems including under-load 
tap changing transformers (ULTC), step voltage regulators 
(SVR), and fixed and switchable capacitors (PC and SC). 
The primary function of these devices is to maintain volt- 
age at the customer service entrance within an acceptable 
range to ensure adequate operation and lifetime of customer 
equipment. In many distribution systems, the operation of 
these devices is governed by local conditions primarily local 
voltage and current sensors. Utilities are increasingly utilizing 
distribution-level Supervisory Control and Data Acquisition 
(SCADA) to operate this equipment, and in some cases, 
this has led to increased levels of centralized control of this 
equipment allowing for coordinated operation. Centralized, 
coordinated control also provides the opportunity to optimize 
the operation to meet utility goals such as minimization of 
losses, reduction of peak apparent power, or extension of 
equipment life. A review of previous work in optimal place- 
ment and sizing of PCs and SCs and coordinated operation 



of SCs and ULTCs can be found series of papers by Baran 
and Wu m], O and Baldick and Wu [3]. 

The relatively slow operation of ULTCs, SVRs, and SCs is 
acceptable in most distribution systems where fluctuations in 
loads and voltage levels are relatively small and significant 
changes in average load occur relatively slowly and in a pre- 
dictable fashion through out the day and year However, this 
situation is due to change as distribution systems are subject 
to higher levels of time-variable distributed renewable gener- 
ation, primarily residential solar photovoltaic (PV) systems. 
High penetrations of time-variable renewable generation will 
pose several new challenges to voltage regulation but may 
also create new opportunities for optimization of distribution 
systems. 

The distributed nature of residential PV generation implies 
that real power is injected at many points along the circuit 
making it no longer possible to obtain reliable estimates of 
the power flows throughout the circuit from a few mea- 
surements of current made at a few discrete locations. In 
addition, real power flows will increasingly be two way and 
in some circumstances, the present deployment of ULTCs, 
SVRs, and SCs, may not be sufficient to ensure adequate 
voltage regulation at the service entrances. ||4| Purthermore, 
the variability of PV generation can occur on a timescale 
much shorter than the present equipment can cope with. Por 
instance, cloud transients can cause ramps in PV generation 
on the order of 15% per second at a particular location 
slowing to perhaps 15% per minute for an entire distribution 
circuit due to its spatial diversity. ||5], ||6| At high levels of 
PV penetration on a distribution circuit, this can result in a 
reversal of real power flow (i.e. change from net generation 
to net consumption) over a period of a few minutes and a 
loss of voltage regulation due to the slow response of existing 
equipment. The existing equipment could respond on shorter 
timescales, however, the increased number of operations 
that would be required to counteract the variability due to 
weather conditions would drastically reduce the lifetime of 
the switches and tap changers ID, ||6l. 

To mitigate many of the issues discussed above, it has 



been proposed that the interconnection standards for inverter- 
based distributed generation be changed in such a way 
to enable the inverters to assist with high speed voltage 
regulation ID, ||6l, Q- To be compliant with present intercon- 
nection standards, fSll PV inverters must not inject or consume 
reactive power or in anyway attempt to regulate voltage. 
However, inverters with this control capabihty already exist 
for off-grid applications and for grid-tied applications when 
the PV system is operating in an islanded mode. Availability 
of hardware is not a fundamental barrier, but excess apparent 
power capacity (above the real power capability of the PV 
system it is connected to) must be built into an inverter 
to allow for reactive power generation and consumption 
while operating near maximum real power. Determining the 
appropriate size of this additional capacity is an important 
outstanding question and depends on yet-to-be-developed 
control schemes that will coordinate the inverters' response to 
changes in voltage and power flow. In addition to effectively 
regulating voltage, a coordinated control scheme should also 
allow for the optimization of the reactive power flows to 
minimize dissipation in the distribution circuit. This paper 
presents a case study of such an algorithm utilizing optimal 
and distributed control of PV-inverter reactive power gener- 
ation. 

The layout of the material in the remainder of this 
manuscript is as follows. Section describes a simplified 
model of an inverter capable of limited reactive power 
generation and consumption. Section |III] describes an opti- 
mization problem where we utilize the inverters' additional 
capacity to minimize losses in a radial distribution circuit 
while respecting the constraints of voltage regulation and the 
inverters' apparent power capacity. In SectionHVl we describe 
the parameters for a model distribution circuit that serves 
as a prototype of a sparsely-loaded rural distribution circuit. 
Section [V] reports the simulation results demonstrating fea- 
sibility of the distributed control and illustrating the quality 
of improvements possible on our prototype circuit. Finally, 
Section [Vl] discusses our conclusions and path forward. 

II. Inverter as a limited regulator of local 

REACTIVE POWER FLOW. 

The present interconnection standard for PV-inverters fHll 
forces all inverters to operate at unity power factor while 
operating in a grid-tied mode, i.e. the inverter must not 
generate or consume reactive power or attempt to regulate 
voltage. Several researchers have proposed that, to allow 
for high penetrations of PV on a distribution circuit, the 
present standard is not workable because the "utility-scale" 
regulation equipment discussed earlier is not sufficient to 
handle two-way power flows or fast enough to mitigate rapid 
cloud transients lH, ||6l, iQ. To be used in this capacity, 
PV inverters must be allowed a new degree of freedom to 
provide voltage regulation. A possible scenario is that a PV 
inverter will be required to have a maximum apparent power 
capability s larger than the maximum power output of its 
PV panel array, maxp^s\ and the excess capability will be 
dispatched by the distribution utility to provide for voltage 
regulation. 



Fig. 1. When s is larger than p^^\ the inverter can supply or consume 
reactive power q^^\ The inverter can dispatch g'^^ quickly (on the cycle- 
to-cycle time scale) providing a mechanism for rapid voltage regulation. As 
the output of the PV panel array p^^^ approaches s, the range of available 
^ decreases to zero. 




A simple model for the relationship between the various 
PV inverter output variables has been described previously. ||7l 
If s is larger than p'-^K the inverter can supply or consume 
reactive power q'^^K The magnitude of q^^^ is bounded by 
■\/s2 — )2 and decreases as the real power output of the 
inverter approaches s. The phasor relationship between the 
inverter operating parameters is shown in Fig. [T]for several 
different levels of PV panel power output, p^^^f . 

Although advanced inverters may have the capability to 
generate and consume reactive power, this output must be 
dispatched in such a way to effectively regulate voltage and 
achieve other utility goals. The algorithms used to perform 
the dispatch could be based on either local conditions or 
be done centrally using circuit-wide information. Because of 
the diversity of distribution utility infrastructure, operating 
areas, and economic models, we develop and compare both 
approaches. 

To deploy centralized algorithms, we envision that the 
communication and control functions of smart-grid technolo- 
gies will be crucial. Enhanced communications capabilities of 
advanced smart-grid systems will allow distribution automa- 
tion systems (DAS) high speed access to voltage amplitude 
and power measurements at many if not all service entrances 
greatly expanding the DAS's grid visibility and situational 
awareness. These measurements, coupled with power flow 
models of the circuit, will allow optimization algorithms 
(developed here and by others) running in the DAS to 
individually dispatch reactive power from each PV inverter 
to ensure that service voltages stay within acceptable bounds 
and, for instance, to minimize losses. With this level of 
system knowledge and centralized approach, algorithms can 
be developed that guarantee optimal solutions. 

However, not all smart-grid schemes allow for high speed 
communication, e.g. power line carrier (PLC) schemes may 
be too slow to update voltage and power measurements fast 



Fig. 2. Diagram and notations for the radial network. Pj and Qj represent 
real and reactive power flowing down tlie circuit from node j, where Pq and 
Qo represent the power flow from the sub-station, pj and qj correspond to 
the flow of power out of the network at the node j, where the respective 
positive [negative] contributions, p^'^^ and g'.'^' [p and g'*'] represent 
consumption [generation] of power at the node, ^he node-local control 
iv) 

parameter qj can be positive or negative but is bounded in absolute value 
as described in Eq. |4] The apparent power capability of the inverter Sj is 
preset to a value comparable to but larger than maxp 



or instantaneous PV generation), while the reactive power 
generated by the PV inverter, qf^ , can be adjusted. However, 
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The rate of energy dissipation (losses) in the distribution 
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is an important global characteristic. Minimizing or at least 
keeping the losses acceptably low is a natural goal for opti- 
mization and control. However, voltage variations along the 
circuit must stay within strict regulation bounds. Measured 
on a per unit basis, the voltage bounds become 
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where normally e w 0.05. 

Combining all of the above, we arrive at the following 
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enough to mitigate the effects of cloud transients. In this 
case, local schemes that only rely upon local measurements, 
perhaps including knowledge of nearest-neighbor nodes, are 
developed. In this article, we develop both an optimal cen- 
tralized and suboptimal but local algorithm and compare their 
performance. The results demonstrate that the performance of 
a simple distributed scheme can approach the performance of 
centralized schemes. 

III. Power Flow. Optimization of Losses and 
Voltage Control. 

To solve for complex power flows in a radial circuit, we 
follow the DistFlow description of H], 121, ||9l- The system 
of AC power flow equations representing the radial circuit 
illustrated in Fig. ^ is Vj — 1, ■ ■ ■ ,n: 



for known 



Eqs. (I1I2I3I4I6I I 

impedances and given configuration of 
py-i ^pyni ^qy-i . The general DistFlow problem is not 
convex and may have multiple solutions. We reduce 
the computational burden by considering the simplified 
DistFlow problem based on the DC approximation, i.e., 
LinDistFlowQ, 0, E 
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where Pj + iQj is the complex power flowing away from 
node i toward node j + Vj is the voltage at node j, Vj + ixj 
is the complex impedance of the link between node j and 
j + and pj+iqj is the complex power extracted at the node 
j. Both Pj and qj are composed of local consumption minus 
local generation due to the PV inverter, i.e. pj = p''^' — p^^'' 
and qj 



q'^p '^q'^f'' ■ Of the four contributions to pj+iqj, p'j 



pY' and, qj'^^ are uncontroUed (i.e. driven by consumer load meters. The real power consumed at each node (p'j'') is 



The formulation is a convex quadratic problem (as the 
quadratic cost function is convex and all the constraints are 
linear) with unique solution that can be computed efficiently. 
We will argue in Section |IV] that this DC-based approxima- 
tion is well justified for our example of a rural distribution 
circuit. 

IV. Description of the prototypical rural 

DISTRIBUTION CIRCUIT 

To demonstrate the technique, we consider a model dis- 
tribution circuit based loosely on one of 24 prototypical 
distribution circuits described in a taxonomy of distribution 
circuits. ifTOl Our model represents a sparsely-loaded rural 
distribution circuit with a nominal line-to-neutral voltage of 
7.2 kV. The line impedance (0.33 + 0.38i)il/ km is constant 
and based on typical conductor types, sizes and spacings. lfTOll 

For this initial study, we consider 100 load nodes separated 
by distances uniformly distributed in the range 200 to 300 



Fig. 3. Energy saved (in percent of the energy lost if all g'^' are set to zero) 
as a function of s. Here, s is reported in kW, and s = 1.1 implies that the 
apparent power capacity of the inverters is 10% higher than the maximum 
real power output of the PV system. The different curves correspond to 
configurations with different levels of PV penetration r. 
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selected from a uniform distribution between and 4 kW, 

(c) 

and the reactive power consumed (qj ) is selected from 

(c) (c) 

uniform distribution between 0.2pj ' and 0.3p^ yielding 
power factors for each load distributed between approxi- 
mately 0.96 to 0.98. At nodes with PV-generation capability, 
the real power generated (pj^^) was always 1 kW reflecting 
a situation where the solar insolation is constant over the 
distribution circuit and all PV systems are the same size. We 
explore different PV penetration levels by randomly assigning 
a variable fraction of the nodes to have PV generation. At 
those nodes, we assign the same apparent power capacity s 
and the reactive power generated (or consumed) is bounded 
by Eq. H) We explore the effects of different amounts of 
inverter excess apparent power capacity by varying s. 

In our model the characteristic values of the nonlinear 
terms cx {Pf + Qf)lyf in Eqs. ( [112] ) are about 10^ times 
smaller compared to the linear terms Pj , Qj , so modeling 
based on the LinDistFlow DC approximation ([8]) produces 
results almost indistinguishable from the exact AC model 

V. Simulations: Results and Discussions 

For a given s and PV-penetration factor r, we considered 
many different realizations of the prototype circuit. Each real- 

(c) (c) 

ization consisted of different % ' , and distances between 
adjacent nodes. The realizations are generated by drawing 
from the distributions described in the previous section. By 
considering many samples for each combination of s and 
r, we observe that the important quaUtative and quantitative 
results described below are sufficiently robust. Therefore, 
in this publication, we present results for a typical sample 
and leave the detailed discussion of statistics for further 
and more formal exploration. Each realization serves as the 
entry point for the DistFlow or LinDistFlow optimization 
algorithm. However, as argued above, the nonlinear terms in 
the AC power flow are quite small and LinDistFlow provides 
an accurate solution. 



Fig. 4. Energy saved (in percentage of the total energy loss observed when 
all g(f ) are set to zero) as a function of PV penetration measured in terms of 
the fraction r of nodes that can inject reactive power. Here, we set s = 1.1. 
The solid-red curve shows results of global optimization. The dashed-green 
curve corresponds to the case of local control where the injected reactive 

power is equal to either the reactive power consumed or to — (p^"')-^ 
if the consumed power is higher than the respective bound on gt^' given in 
Eq.|4] 
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In our preliminary study of this problem, we investigate the 
reduction in losses first as a function of s for several different 
values of r. Second we investigate the reduction in loss as a 
function of r for a single value of s. We compare this result 
to that obtained via a simple distributed control approach. 
Finally, we present the improvement in power quality (i.e. 
voltage regulation) obtained by using a centralized algorithm 
to dispatch gj^"* from the nodes with PV capability. 

Figure [3] shows that the energy savings increase monotoni- 
cally with the inverter apparent power capacity s. For a wide 
range of PV penetration, the majority of the energy savings 
occurs by s = 1.1, i.e. at an apparent power capacity of only 
10% higher than maxpj^''. Additional savings are possible 
beyond s ~ 1.1, however, higher values of s will increase 
the cost of PV systems. Determining the optimal value of s 
must also include economic factors and is beyond the scope 
of this work. The energy savings and the value of s at which 
it saturates will depend on the reactive power consumed by 
the loads, however, we have not explored this dependence 
in this preliminary study. Also, we have only considered the 
case when all the pj^-* are at their maximum value. We note 
that the saturation observed in Fig. [3] will occur faster in 
more realistic models of renewable generation where the PV 
systems will not always operate at their maximum capacity. 

Figure H] shows how the energy savings depend on PV 
penetration for s = 1.1, i.e. a vertical slice through Fig [3] at 
s = 1.1. First we note that the energy savings approaches 
20% at very high PV penetrations suggesting that, if reactive 
power dispatch can be achieved, high PV penetration is in 
fact beneficial. Again, we note that the quantitative results 
presented in this manuscript depend on our assumptions of 
the reactive power consumed by the loads. However, we do 
not expect the qualitative behavior to change as the load 
assumptions are changed. Figure |4] also compares a global 
optimization (red) with a local control scheme Q requiring 



Fig. 5. Voltage variation and qj® for a particular circuit realization with 
T = 90% and s = 2.0. The voltage drop (normalized with respect to initial 
voltage Vq) is less if one is allowed to inject reactive power into the circuit. 
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no communication at all. In this particular local scheme, g^^^ 
is set to minimize the local reactive power consumption qj 
subject to the constraints of Eq.|4] For the prototypical circuit 
and load assumptions used in this preliminary study, this 
naive control strategy surprisingly achieves about 95% of the 
maximum possible savings. The efficacy of this local scheme 
should be explored for many different circuit configurations 
and loads. Note that both strategies should be compared with 
the "do nothing", qj^'' = case shown in dashed black in 

Fig. a 

Finally, we demonstrate the effect on power quality (i.e. 
voltage regulation) of dispatching reactive power from a 
relatively high penetration of PV inverters. Figure |5] shows 
the voltage profile along the circuit with (green) and without 
(red) reactive power dispatch. Here, the reactive power is 
dispatched using the global scheme described above. In 
addition to significant reduction in losses, the controlled 
dispatch of reactive power increases the power quality by 
significantly reducing the voltage drop along the circuit. 
Although not shown in the Figure, we note that the global 
control of reactive power performs better than the local 
control because it allows for better suppression of reactive 
power flow Qj over the whole line, not just qj. Figure |5] also 
shows the individual q'^^^ at each node along the line (vertical 
blue dashed lines, each standing for a generator with some 
lines missing because r < 1.0). 

VI. Conclusions and Path Forward 

Our preliminary study is meant to elucidate the qualitative 
behavior of a distribution circuit with various levels of PV 
system penetration, excess inverter apparent power capacity, 
and global versus local inverter-reactive power dispatch. We 
find that power dissipation is reduced and power quality 
is increased as both PV penetration and excess apparent 
power capacity are increased. For the prototypical circuit and 
loads considered in this study, we find that the reduction in 
power dissipation plateaus at a relatively low value of excess 
apparent power and that a local control scheme performs 
nearly as well as a global solution. 



Perhaps more importantly, this study helps to formulate 
future research directions to improve the optimization and 
control algorithms. We envision the following future modifi- 
cations to the basic scheme discussed in the manuscript: 

• Different circuit configurations and load profiles must be 
explored to determine if the qualitative and quantitative 
results apply to a wide variety of cases. The taxonomy 
of distribution circuits described in ifTOI provides a guide 
to selecting additional case studies. 

• The global optimality of the algorithm can be traded for 
local optimization. The naive local scheme discussed 
in the manuscript, which requires making only local 
decisions, is already significantly better than no compen- 
sation at all. However, there are all other intermediate 
strategies which may close the gap between the fully 
local (simple but strictly suboptimal) and fully global 
(optimal but complex and requiring significant commu- 
nication and coordination) strategies. We anticipate that 
a linear-scaling algorithm capable of efficiently solving 
the single-branch DistFlow optimization problem may 
be available. We envision this type of optimal algorithm 
to have a form of dynamical programming with 
the global optimization replaced by a sequence of local 
optimizations advanced sequentially from the end of the 
circuit towards the entry-point (and/or in the reversed 
direction). These techniques could also be extended 
beyond the single-branch circuit to the case of a circuit 
with multiple branches, provided that the circuit remains 
a tree (no loops). 

• Accounting for loops in the graph constitutes an inter- 
esting algorithmic challenge that may have applications 
for highly meshed distribution circuits typical of urban 
systems. In this case, solving the DistFlow equations 
will likely be of exponential complexity (in the number 
of nodes), however the LinDistFlow equations are still 
expected to be polynomially-tractable, i.e. reducible 
to a linear program. Normally, linear programming is 
not distributed, however in some special cases, e.g. 
minimum-cost network flow lfT2l . it allows a distributed 
implementation via Belief Propagation (BP) algorithms. 
Developing a distributed BP-algorithm for LinDistFlow 
optimization over a loopy graph is thus another chal- 
lenging task for future exploration. The distributed na- 
ture of the algorithm will naturally enable the ability 
to carry out the global optimization via a sequence of 
local computations and communications limited only to 
neighboring loads on the circuit. 

• Following the general arguments of ifTSl . adding some 
switching capabilities can have a beneficial impact on 
distribution circuit losses for tree-like circuits and po- 
tentially for loopy circuits. Analysis of both algorithmic 
and phase transition aspects of switching for such a 
loopy circuit can be carried out in the spirit of llT4l . 

na. 
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